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ABSTRACT 

A powerful experimental approach for investigating computation in networks of biological neurons is 
the use of cultured dissociated cortical cells grown into networks on a multi-electrode array. Such prepa- 
rations allow investigation of network development, activity, plasticity, responses to stimuli, and the 
effects of pharmacological agents. They also exhibit whole-culture pathological bursting; understanding 
the mechanisms that underlie this could allow creation of more useful cell cultures and possibly have 
medical applications HI 13. 



This paper presents preliminary results of a computational study of the interplay of individual neuron 
activity, cell culture development, and the network behavior. We investigate whether bursting can occur 
in an initially unconnected "network" that develops connections according to an experimentally-verified 
model of cell culture connectivity growth. 



Neuron Model An integrate-and-fire neuron model with dynamical synapses that exhibit activity- 
dependent facilitation and depression was used |l3l|4l|5l. It includes synaptic, constant, and noise currents, 
with a reset of its membrane voltage to a fixed value upon exceeding threshold and generating a spike 
and a fixed absolute refractory period thereafter. 

The synapse has four state variables: three that govern the fraction of synaptic resources in particular 
states — X (recovered state), y (active state), and z. (inactive state) — and one, u, that represents synaptic 
efficiency (see the more complete paper* for detailed equations). See 10 for parameter values used. 



Network Model Simulations were conducted by constructing networks with model neurons on a rect- 
angular grid. Connectivity was determined by incorporating a model of cortical cell culture connectivity 
development m that model's neurite outgrowth as a radius of connectivity that changes at a rate inversely 
proportional to a sigmoidal function of cell firing rate: 

-^=PO(F,) (1) 

"'^''^'- l+e.p(('-.0//i) ® 

where /?, is the radius of connectivity of neuron /, Fi is neuron /'s normalized firing rate, p is a rate 
constant, £ is a constant that sets the "null point" for outgrowth (the normalized target firing rate that 
causes no outgrowth or retraction), and jS determines the slope of G(-). One of the parameters varied in 
these simulations was e. 

Synaptic strengths Wi were computed for all pairs of neurons that had overlapping connectivity regions 
as the area of their circles' overlap. The bulk of the neurons in the network were excitatory; a small 



* A more complete version of this paper is available at http://faculty.washingon.edu/stiber/Public/NC07.pdfl 
^ To whom correspondence should be sent. 
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Figure 1: Firing rates normalized relative to target rate (A); burstiness index (B). Both computed for 
25,000-30,000sec and plotted versus the two simulation parameters. 



number were chosen to be inhibitory. Similarly, most neurons were not spontaneously active, but a few 
had lowered firing thresholds to produce spontaneous firing at a rate of between 0.02 and 6 spikes/sec. 
To produce more consistent simulation results, a set of standardized layouts was chosen to maximize 
spacing among inhibitory and spontaneously active cells and reduce edge effects. The fraction of cells 
that were inhibitory was the other simulation parameter varied. 

Computer Implementation We used CSIM (A Neural Circuit SIMulator) version 1.1 for the simula- 
tions. The original code was pared down to a small core that was linked to a stand-alone C-i-i- program to 
run on Linux, Windows, and Macintosh computers. Generally speaking, each simulation took between 
two and 20 hours on computers with 2-3GHz microprocessors. 

Simulations consisted of networks of 100 neurons in a 10x10 arrangement simulated for 30,000-60,000 
seconds. While the rate of neurite outgrowth was greatly sped up compared to the living preparation, 
numerical investigation indicated that this did not introduce instability in the simulation. 

Analysis Methods To examine global behavior, average neuron firing rate and burstiness index (BI) O 
were calculated and plotted versus the two simulation parameters. BI was computed by first calculating 
a spike count vs. time histogram for the entire network during the last 5,000sec of the simulation. 
The fraction, /15, of the total number of spikes contained by the 15% most populous bins was then 
normalized to produce the burstiness index, BI, as BI = (/15 — 0.15)/0.85. Detailed examination of 
single simulations involved plots of neurons' connectivity radii and firing rate versus time. 

Results An initial set of 50 simulations was performed with target rates in the range 0.1-1.9 (inclusive, 
in 10 steps) and fraction of excitatory neurons 0.9-0.98 (inclusive, in five steps). Figure [TjA) shows 
normalized firing rates for the final 5,000sec of each simulation. Only the simulations with higher target 
rates and fewer inhibitory neurons showed the great increase in firing rate that might be associated with 
bursting. Longer (60,000sec) simulations exhibited some bursting type of behaviors with the fraction of 
bursting as low as 0.9 for the higher target rates. Figure [TjB) indicates that the fraction of excitatory 
neurons in this range has only a modest effect on burstiness index (for 30,000sec simulations), and that 
moderate target rates produce the highest BI values. This apparent conflict with the previous observation 
of bursting at higher target rates can be explained by examining the detailed behavior of individual 
simulations. 

Figure |2] shows detailed information for simulations with three sets of parameters: (target rate, fraction 
excitatory cells) = (0.1, 0.9), (0.9, 0.94), and (1.9, 0.98). These include the parameter extremes and a 
central value and both bursting and non-bursting activity. In these cases, 60,000sec simulations were 
performed. Note that some of the simulations that weren't bursting at 30,000sec were bursting shortly 
thereafter, as evidenced by the (0.9, 0.94) one. 
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Figure 2: Detailed simulation results for simulations with parameters (target rate, fraction excitatory 
cells) of (0.1, 0.9) (left), (0.9, 0.94) (middle), and (1.9, 0.98) (right). Data for edge neurons are green, 
non-edge (and non-inhibitory, non-spontaneously active) neurons are black, inhibitory neurons are red, 
and spontaneously active neurons are blue. 
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Figure 3: Effects of neurite outgrowth on bursting for simulations with parameters (target rate, fraction 
excitatory cells) of (0.9, 0.94) (left, middle) and (1.9, 0.98) (right). Growth was "frozen" either during 
(left, right) or between (middle) bursts. 



Nevertheless, this confirms that the low BI values for low target rates correspond to non-bursting behav- 
iors (and that the connectivity radii had not stabilized for cells that were not inhibitory or spontaneously 
active). Low BI values for high target rates were a possible result of the very broad or frequent bursts. 
For the bursting behaviors, connectivity radii have stabilized, excepting small variations during bursting. 
In all simulations, connectivity radii for edge neurons are larger than others, inhibitory neurons had mod- 
erate radii, while spontaneously active neurons had a wide range of different connectivities, likely due to 
the variability in their firing thresholds. 

In either bursting or non-bursting behaviors, spontaneously active neurons tended to be the most active. 
This is not surprising, as their lowered thresholds would make them more excitable. The next higher 
firing rates belonged to the inhibitory cells, then non-edge cells, then edge cells. 

In the results in figure |2] it seems possible that the mechanism for burst initiation and termination is the 
variation in connectivity. Figure [3] presents simulations in which connectivity was fixed either during 
bursts (left, middle) or in between bursts (right). For the lower target rate/lower fraction of excitatory 
cells simulation, it does indeed seem that bursting is controlled by connectivity. However, for a higher 
rate/higher fraction of excitatory cells, bursting can continue even in the absence of variation in connec- 
tivity (this simulation also produced bursts when growth was stopped during a burst). 

Discussion Bursting occurred with sufficiently small inhibition and high target firing rate. One might 
expect the latter to produce greater connectivity for every neuron, which in turn would be the mecha- 
nism for whole-culture bursting. However, in the low-target-rate, non-bursting simulations, such as fig- 
ure [2jleft), non-inhibitory, non-spontaneously active cells grow large connectivity radii. Spontaneously 
active cells, on the other hand, tend to have large connectivity radii in the bursting simulations. Presum- 
ably, lowering these cells autonomous firing rate would result in bursting at lower target rates. 
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In previous investigations of bursting with randomly connected networks f4], the model synapses' de- 
pression and facilitation were, neglecting the influence of noise, the mechanism underlying burst ini- 
tiation and termination. Our preliminary results indicate that this is possibly the case under certain 
circumstances, but not all. For some regions of parameter space, it may be the case that the mechanism 
is a hysteresis effect involving changing connectivity radii. There are a number of possible reasons for 
this difference: 

• A number of parameters were set arbitrarily or not fully explored. These include fraction of ex- 
citatory cells, scaling of synaptic weights from area of connectivity overlap, and no differential 
scaling based on type of synapse (i.e., inhibitory vs. excitatory). 

• In the current simulations, only spontaneously active cells had any parameter variability; all other 
cells of a given type (inhibitory or excitatory) had identical parameters. 

• The overriding issue here is likely the small network size. Edge effects were great (edge neu- 
rons' connectivity radii were always the greatest of all neurons and 36/100 of the cells were edge 
neurons) and there were small numbers of inhibitory and spontaneously active cells. The final net- 
works were almost completely connected. Increasing network size to, say, 100x100, could have 
little impact on final connectivity radii but with each neuron having connections to less than 10% 
of the network. 

Increasing network size will have computational consequences that must be addressed: in its cur- 
rent form, a 60,000sec simulation of a 100x100 network would take at least 2,000 hours (83 days). 

There are also fundamental differences between the connectivity patterns generated by this model (per- 
haps most similar to radial basis functions) and many other models of cortex or recurrent networks (in 
which either network topology is irrelevant or a power law-type distribution is used that produces mostly 
local connections with a few long-range ones). It will be instructive to investigate the detailed correlation 
structure of inter- and intra-burst neuron firing. 

Keywords: cortical cultures, network development, bursting. 
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